Trojan Horse as an indirect technique in nuclear astrophysics. Resonance reactions. 
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The Trojan Horse method is a powerful indirect technique that provides information to determine 
astrophysical factors for binary rearrangement processes x + A — > b + B at astrophysically relevant 
energies by measuring the cross section for the Trojan Horse reaction a + A — *• y + b + B in quasi- free 
kinematics. We present the theory of the Trojan Horse method for resonant binary subreactions 



■ based on the half-off-energy-shell R matrix approach which takes into account the off-energy-shell 

' effects and initial and final state interactions. 

bJO! PACS numbers: 26.20.+f, 24.50.+g, 25.70.Ef, 25.70.Hi 
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I. INTRODUCTION 



The presence of the Coulomb barrier for colliding charged nuclei makes nuclear reaction cross sections at astrophys- 
ical energies so small that their direct measurement in the laboratory is very difficult, or even impossible. 

Consequently indirect techniques often are used to determine these cross sections. The Trojan Horse (TH) method 
is a powerful indirect technique which allows one to determine the astrophysical factor for rearrangement reactions. 
The TH method, first suggested by Baur [l[ , involves obtaining the cross section of the binary x + A —> b + B process 
at astrophysical energies by measuring the two-body to three-body (2 — > 3) TH process, a + A^y + b + B, in the 
quasi- free (QF) kinematics regime, where the "Trojan Horse" particle, a = (xy), is accelerated at energies above the 
Coulomb barrier. After penetrating through the Coulomb barrier, nucleus a undergoes breakup leaving particle x to 
interact with target A while projectile y flies away. From the measured a + A — > y + b + B cross section, the energy 
00 . dependence of the binary subprocess, x + A — > b + B, is determined. 

The main advantage of the TH method is that the extracted cross section of the binary subprocess does not contain 
the Coulomb barrier factor. Consequently the TH cross section can be used to determine the energy dependence of 
the astrophysical factor, S(E), of the binary process, x + A — ► b + B, down to zero relative kinetic energy of the 
particles x and A without distortion due to electron screening [1,[3[. The absolute value of S(E) must be found by 
normalization to direct measurements at higher energies. At low energies where electron screening becomes important, 
comparison of the astrophysical factor determined from the TH method to the direct result provides a determination 
of the screening potential. 

Even though the TH method has been applied successfully to many direct and resonant processes (see [|| and 
references therein), there are still reservations about the reliability of the method due to two potential modifications 
of the yield from off-shell effects and initial and final state interactions in the TH 2^3 reaction. Here we will address 
the theory of the TH method for resonant binary reactions x + A — > b + B. 



II. TROJAN HORSE 

The TH reaction is a many-body process (at least four-body) and its strict analysis requires many-body techniques. 
However some important features of the TH method can be addressed in a simple model. Let us consider the TH 
process assuming that nuclei y, x and B are constituent particles, i. e. we neglect their internal degrees of freedom. 
For simplicity, we disregard the spins of the particles. The TH reaction amplitude is given in the post form by 

M(P,Ka) =< xfci 4 _) |AVl*l +) > ■ (1) 

Here, is the exact a + A scattering wave function, $^T^ is the wave function of the system F = b + B = x + A, 
Xk p( r u) ^ s ^ ne distorted wave of the system y + F, ipi is the bound state wave function of nucleus i, r,y and ky are 
the relative coordinate and relative momentum of nuclei i and j, P — {k y F,^bB} is the six-dimesional momentum 
describing the three-body system y, b and B in the final system, AV v f — V v f — UyF, V v f = V y b + V v b — V yx + V v a 
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is the interaction potential of y and the system F and U yF is their optical potential. The surface approximation 
suggested in [? ] was the first serious attempt to address the theory of the TH method. The surface approximation 
assumes that the TH reaction amplitude has contributions from the external region where the interaction between the 
fragments b and B (x and A) can be neglected and the wave function <I>^ ' can be replaced by its leading asymptotic 
form 



*<+> « Vb [e* k — +F bB u{+l(r bB )] + J^M bB ^ xA ^JL- ««(r , , ,. 



where $p = ^tafF) an< ^ ^L^bCF) = -k* B (F)' ( r y) * s ^ ne outgoing spherical wave, F bB is the b + B elastic 
scattering amplitude, M bB ^ xA is the b+B — > x+A reaction amplitude inverse to the binary reaction x+A — > b+B and 
My is the relative velocity of nuclei i and j. The expression for the TH reaction amplitude in the surface approximation 
is given by 

M(P,k aA ) ~ M 6B ^ A < X { kv i ^Aui^(r xA )\AV yF \(p a ( f Ax'C ) A (r aA ) >, (3) 

where the exact initial scattering wave function is replaced by ip a ip A Xktl ( r a^) an< ^ xIa ^ s tne distorted wave 
describing the scattering of the nuclei a and A in the initial state of the TH reaction. For simplicity we don't take 
into account here the Coulomb interactions. However, in the case of the resonant binary reaction x + A — > b + B the 
dominant contribution comes from the nuclear interior where both channels x + A and b + B are coupled and where 
the asymptotic approximation for cannot be applied [l2l|. 

In this work we will address the theory of the TH method for the resonant binary subprocesses x + A — > b + B 
which explicitly takes into account the off-shell character of x. Eq. ([T|) can be used as a starting point to derive the 
expression for the TH reaction amplitude. We assume that the resonant reaction x + A — > b + B proceeds through the 
formation of the intermediate compound state <&i, i. e. we neglect the direct coupling between the initial x + A and 
final b + B channels, which contributes dominantly to direct reactions but gives negligible contribution to resonant 
ones. An important step in deriving the resonant contribution to the TH reaction matrix element is the spectral 
decomposition for the wave function &p ' given by Eq. (3.8.1) 0- It leads to the shell-model based resonant R 
matrix representation for ' which is similar to the level decomposition for the wave function in the internal region 
in the R matrix approach: 

N 

E K 6B (^)[D-V^r- (4) 

V,T=1 

Here N is the number of the levels included, E bB is the relative kinetic energy of nuclei b and B, $ r is the bound 
state wave function describing the compound system F excited to the level r. D VT is similar to the level matrix in 
the R matrix theory and is given by Eq. (4.2.20b) 0). Finally, 

V* B {E bB ) =< X ( bB ] w|AV 6 b|*„ > (5) 

is the resonant form factor for the decay of the resonance F v described by the compound state into the channel 
b + B. The partial resonance width is given by 

t v {E bB ) =2Tr\V» B (E bB )\ 2 . (6) 

Then the TH reaction amplitude is 

N 

M^(P,k aA )^ Yl V^iEhB)^- 1 }^ M r (k yF ,k aA ), (7) 

where M T (k y p, k aA ) is the exact amplitude for the direct transfer reaction a + A — ► y + F r populating the compound 
state F T of the system F = x + A = b + B: 

M T (k yF ,k aA ) =< x yF ] $ T \AV vF \*l +) > . (8) 

The direct transfer reaction is very well described by the DWBA amplitude, i. e. for the practical analysis we can 
approximate ~ <p a tp A X^a - Correspondingly, M T (k yF ,k aA ) can be replaced by 

M? W (k yFl k aA ) =< Xy'J QrlAVyFlcpaifiAX^ > ■ (9) 
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Correspondingly for the TH reaction amplitude we get from Eq. 

N 

M( R \P,k aA ) a VZ B {E bB )[T>-X T M? w {k yF ^ aA ). (10) 

i/,r — 1 

The DWBA amplitude takes into account the rescattering of nuclei a and A in the initial state of the TH reaction 
and enters as a form factor into the TH resonant reaction amplitude reflecting the off-energy shell character of the 
transferred particle x. Since in the TH method the astrophysical factor determined from the TH method is normalized 
to the on-energy-shell (OES) S factor, the replacement of the exact transfer amplitude by the DWBA one, as we will 
see, practically does not affect the final result. 



A. Single resonance 



The triple differential cross section for the TH process a + A — > y + b + B proceeding through an isolated resonance 
F T is given by 

d 3 a _ . T bB{T) (E bB )\M T DW (k yF ,k aA )\ 2 

A 3 — — W7fT7-\ ■ \ LL ) 



dE bB dn kbB dn KF {ExA -e Rt )z + 

Here, A3 is the kinematical factor, T bB ^(E bB ) is the observable resonance partial width in the channel b + B, 
F T (E xA ) is the total observable width of the resonance F T . Note that all functions T(E) are related to T(E) as 
T(E) = T(E)/(1 — ( d ^ B T )e=e Rt )i where A Tr is the t level shift. Also Er t is the resonance energy of the resonance 
F T in the channel x + A. Thus the TH triple differential cross section, in contrast to the OES single-level resonance 
cross section, contains the generalized form factor \M^ W (k y p, \t aA )\ 2 rather then the entry channel partial resonance 
width T xA (t){Exa) of the binary process x + A — > b + B. A simple renormalization of the TH triple differential cross 
section allows us to single out the OES astrophysical factor for the resonant binary subprocess x + A — > b + B: 

S(E xA ) = NF(E xA ) = _!L_ ^(E^T^E^) 

dE bB dr> kbB dn UyF 2fi xA (e xA _ Ert )2 + r T (fixA) 

where the normalization factor NF(E xA ) is given by 

AT Ft TP A - - 1 TP . 2*VxA ^xA(r)(E xA ) 

N nE xA )- WA - 3 E xA e lM o W{kyFMr (13) 

Note that the DWBA amplitude M DW (ky F , Ka) remains practically constant on the interval of a few hundreds keV. 
Eq. p^l> explaines and justifies the phenomenological procedure used before successfully in the TH analysis (see [3| 
and references therein). The renormalization factor can be rewritten as 

NF{E xA )=e 2 * [v * A - V °<^ ] T ^§^ NF(E Rl ), (14) 

where ^ xA (t)(E xA )/^ xA ( t )(Er 1 ) — P(E xA ) is the barrier penetration factor appearing in the R matrix theory. The 
factor NF(Eji 1 ) can be found phenomenologically by comparing the experimental TH triple differential cross section 
with the available OES experimental astrophysical factor at resonance energy. This phenomenological normalization 
leads to the intermediate astrophysical factor 

S'(E A ) = 11 e 2 "^I^A=^ r bB{T) {E bB )r xA{T) (E Rl ) 
The final astrophysical factor can be derived by multiplying S'{E xA ) by the energy-dependent factor in Eq. (|14|) 

2 7T A — A I ] 

e e *a= e r 1 r xA r T \(E xA )/T xA r T )(E Rl ). Thus normalization of the triple TH differential cross section to the 

experimental astrophysical factor at resonance energy achieved by multiplying Eq. (fTTj) by the factor NF(E xA ) plays 
a very special role in the TH method. 
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B. Two interfering resonances 



For two interfering resonances we need to consider the two-level, two channel case. This requires the half-off-cnergy- 
shell (HOES) R matrix formalism. Here we address this formalism for a simple case when the distances between two 
resonances are significantly larger then their total widths. Then the OES reaction amplitude in the R matrix formalism 
is given by the sum of the amplitude of each resonances (see Eq. (XII, 5. 15) 0]). The corresponding expression for the 
HOES reaction amplitude can be obtained by the replacement of the resonance partial widths in the entry channel of 
the binary reaction x + A ^ b + B by the corresponding generalized form factors M® w (k yF , k a ), t = 1,2. Thus the 
triple TH cross section in the presence of two interfering resonances in the subsystem F = x + A = b + B is given by 

d 2 a |V Tl^ {T) {E bB )M? w {), yF ^ aA ) 

dE bB dn kyF dn kbB 3 1 t ^ 2 ExA -e Rt+1 J^f-i 1 ' ( j 

We assume that E Rl < Er 2 . The goal of the THM is to determine the energy dependence of the astrophysical factor 
at the astrophysically relevant energies. The ratio Al^ w — M® w (k yF , k aA )/M[ >w (k yF , k a A) is practically constant 
in the interval of a few hundred keV, E x a < E Rl . Normalizing the TH cross section to the OES S factor at E = E Rl , 
where the contribution from the second resonance can be neglected, gives the astrophysical factor determined from 
the TH reaction 

, TH(F , iggj r (F )U rffg)^) ^ B \ 2) { EbB )M^ 

o {E xA ) = — 1 xA(i)\E x a)\[ twp s H tTFTtJ ■ \ n ) 

W n ExA -E Rl +i^^l ExA ~E R2 +i^A± 11 

This astrophysical factor is to be compared with the OES astrophysical factor determined from direct measurements 

ne^^ A |r ^ {1) {E bB ) r 1 ^ 2 (E bB ) 1(xA)21 
o{E xA ) — — l xA (i){E x a) I ^77? — H pTT? — H • (18) 



Here, 7( x ^)2i = l{xA)2/l(xA)i = T xA (2) < y E xA)/r 1 xA{1) {E xA ) and 7(zA)t 



is the reduced width for the r-th resonance in 



the channel x + A. Each amplitude M® w (k yF , k a ) is complex, but the ratio M^ w may have a small imaginary part. 
The normalization of the TH S factor to the OES one at resonance energy plays a crucial role in the TH method. 
After such a normalization, we need to know only the ratio of the DWBA amplitudes to calculate S TH (E xA ). 



1. Plane wave approximation 

Ratio M^ w can be approximated by the ratio of the corresponding amplitudes calculated in a plane wave ap- 
proximation, because a simple plane wave approximation gives similar angular and energy dependence as the DWBA 
but fails to reproduce the absolute value. It explains why a simple plane wave approximation works well in the TH 
analysis Note that in the plane wave approximation Mf w (k yF , k a A) is replaced by 

M°(k yF ,k aA ) =< e tU ^- r ^ip y <I> T \V y A + V X A\<Pa<PAe l ^ A - T « A > . (19) 

Note that the post and prior forms are equivalent but the post form is more convenient for our purpose. In the QF 
kinematics for sufficiently high momentum of the the projectile A it will interact dominantly with the fragment x 
while the contribution of the term with V v a is minimized. That is why in what follows we neglect the term containing 
V y A- Then the transfer reaction amplitude in the plane wave approximation takes the form 

M°(k yF ,k aA ) «< e * k «™ Jfr | < V xA > xA | I yx e ik ^ A >, (20) 

where I xA =< ip A tp x \Q T > is the overlap function of the wave function of the resonance state F T and the bound state 
wave functions of A and x, Iy X =< tp y (p x \ ip a > is the overlap function of the bound state wave functions of nuclei 
a, x and y, and ip x , < V xA >=< ip A ip x \V xA \ip x ip A >■ The plane wave amplitude M®(k yF , k aA ) can be written in a 
factorized form 

M T (ky F ,k aA ) = [W^ikA -^k F )Yl« x {ky ~ ^k a ). (21) 

m F " m a 

Here, Iy X (p yx ) is the Fourier transform of the overlap function Iy X (r yx ) and 

W?Z(Ka) =< e l ^ A ^ A \ < V xA > xA {y xA )\I xA {v xA ) > (22) 
is the vertex form factor for x + A — > F. T . Then Eq. (|16p for the TH triple differential cross section takes the form 
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FIG. 1: Comparison of the calculated astrophysical factor S TH (E) for 15 N(p, a) 12 C (solid line), where E = E x a, with the 
direct data [lEEl- 
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r^ 2 (r) (^)[^(k A -^k F )] 



dE bB df2 k „ dO k 



= A 3 |/« (k,-^k a ) 



T=l,2 



Ft, 



E x a - E Rt + i 



(23) 



Now we can get the HOES cross section for the binary subprocess x 
section 



da 



dfl r 



HOES 



d 2 a 



dE bB dfl kvF dft 



A — » b + B from the triple differential cross 



(24) 



a 3 \i^ vx )r 



where Tp yx = k y — — k a . Eq. (|24p explains and justifies the procedure used in IA [J] to connect the triple and binary 
TH cross sections. Note that in a strict approach the triple differential cross section is expressed in terms of the 
overlap function J° x rather then the two-body bound state wave function ip a . 



Note that I yx and ip a are related by 



r 



cl/2 



<Pa, 



(25) 



1 /2 

where S yx is the spectroscopic factor. The binary reaction HOES cross section is only intermediate result. The final 
goal is the TH astrophysical factor which can be determined by normalization of the triple differential cross section to 
the OES astrophysical factor in the first resonance peak and is given by Eq. (fTT)) . In the plane wave approximation 
M^ w is replaced by 



[w x F X (k A -^k F ) 

[ W xa(^a - k F ) 



(26) 



If ~ 7(xA)2ii the astrophysical factor S th (E x a) reproduces the OES S factor S{E x a) at energies E x a < 
Eb^. In Fig. [T]the astrophysical factor S th {E x a) for 15 N(p, a) 12 C calculated using Eq. (fT7|) for the TH reaction 
15 N(d,nd) 12 C is compared with the experimental S(E x a) obtained from direct measurements. There are two 1~ 
interfering resonances at E^ = 312 keV and Er 2 = 962 keV. The best fit has been achieved for T x a(i) = = 1-1 



keV, r 6B (i) = r Q( i) = 93.4 keV, r xA(2 ) = r p(2) - 95.31 keV and T bB 



(2) = r Q(2) = 45 keV. To find M% x we used Eq. 

(122p in which the overlap function 1^ is approximated by a single-particle 15 N— p wave function in the Woods-Saxon 
potential calculated in the internal region by a procedure similar to that used in R-matrix method to calculate the 
level eigenfunctions. We find that M 21 « 1.13 while r f( x A)2i — 1-1 ± 0.1. It explains why the calculated S th (E x a) 
shown in Fig. [T]is in an excellent agreement with the direct data. 

We presented the expression for the resonant S factor determined from the TH reaction taking into account the off- 
energy-shell effects within the HOES R matrix formalism and justified a simple plane wave approximation. Validating 
this makes it clear why the TH method is such a powerful indirect technique for nuclear astrophysics. 
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